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Abstract. Theoretical predictions for lateral distribution 
function of electrons in extensive air showers based on scal- 
ing formalism are presented. Our results are tested by com- 
parison with AGASA experimental data taking into account 
the contribution of low-energy muons and simultaneously the 
effect of scintillation detectors response, according to recent 
simulation results. The possibility for ultrahigh energy cos- 
mic ray primary composition deduction from the shape of 
LDF is discussed in detail. 



1 Introduction 

One of the basic characteristics of extensive air showers 
(EAS) of superhigh energies, that can be measured quite ac- 
curately by large ground-based air shower arrays is charged 
particle local density at various distances from the core lo- 
cation. Reliable theoretical predictions on the lateral distri- 
bution function (LDF) of different EAS components in wide 
radial distance range are therefore of great importance in cos- 
mic ray research for both physical interpretation of existing 
experimental data and new experiments design studies. Un- 
fortunately, the problem of fast, proper and at the same time 
adequate to numerous specific experimental conditions cal- 
culations of LDF at large distances from the shower core 
(r > 1 km) is not solved yet. 

Direct Monte-Carlo technique is of course preferable but 
takes unreasonably much computation time. One of widely 
used approaches implements an analytical description of 
electromagnetic subshowers, usually based on different well 
known modifications of NKG formula (L. G. Dedenko et al., 
1975; A. A. Lagutin et al., 1979; A. V. Plyasheshnikov et al., 
1988) obtained for distances up to several hundreds meters 
from the core position. In this case formal extrapolation to 
r ^ 1 km and farther is a source of mistakes. The weak 
point of hybrid methods is the necessity of complete recal- 
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culation of prerecorded libraries of subshowers in case of 
any changes in input assumptions. Established thinning al- 
gorithms seem very promising, but still not sufficiently ef- 
fective in case of very large radial distances, because small 
fraction of distant particles is tracked with very large weights 
(A. M. Hillas, 1997). There are evidences [see, for example, 
A. Filipcic et al. (2000)], that it is possible to improve sta- 
tistical significance of thinning at large distances by artificial 
limiting the weight, which particle can obtain, but such tech- 
nics are still in the stage of development. 

On the other hand, present situation concerning com- 
parisons of lateral distributions measured by scintillation 
counters at Akeno and Yakutsk with each other as well 
as with theoretical predictions at ultrahigh energies (_Eo > 
lO'^^ eV) remains mostly unsatisfactory. Though in sev- 
eral recent publications [see, for example, A. V. Glushkov 
et al. (1999); M.Nagano et al. (2000)] reasonable agree- 
ment between the shapes of experimental and calculational 
charged particle LDFs was achieved, the absolute values 
demonstrate significant discrepancies. As a consequence 
shifting theoretical densities vertically with a factor ~ 1.5 
is now widespread technique utilized during comparisons of 
lateral distributions. Furthermore latest theoretical studies 
(T.Kutter, 1998; A. A. Lagutin et al., 1999; A. 1. Goncharov 
et al., 2000; M. Nagano et al., 2000) indicate strong influence 
of the effect of scintillation detector response on the shape of 
LDF. It is also worth to mention peculiarities of lateral distri- 
butions of charged particles and muons al E^) > b ■ IQ^^ eV 
founded by Yakutsk group, that have not been confirmed 
by other experiments at this moment and can not be ex- 
plained by simulations without exotic physical assumptions 
(A. V. Glushkov et al., 1997, 1999). 

In this paper we check the validity of one-parametric scal- 
ing representation of lateral distribution of electrons estab- 
lished in our earlier works by comparison with AGASA ex- 
perimental data. In order to perform adequate comparisons 
with experiment recent CORSIKA simulation results about 
contribution of low-energy muons and energy deposition in 
scintillation counters are used. We also examine here the 
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possibiUty of utihzing the experimental data about the shape 
of charged particle LDF for cosmic ray primary composition 
deduction. 



2 Scaling formalism for LDF in cascade showers 

According to scaling formalism (A. A.Lagutin et al., 1997, 
1998, 1999), the lateral distribution of electrons in both 
gamma- and hadron-induced cascade showers can be repre- 
sented in a form: 



Pe{r;Eo,t) 



Ne{Eo,t) 



(1) 



Here pe{i"; E^J) - electron density at radial distance r from 
the core in showers of primary energy _Eo at depth t from pri- 
mary particle injection point, A^e - total number of electrons 
at observation level, -Rm.s. - mean square radius (the second 
moment of normahzed electron LDF), which is defined as 



-Rm.s.(-Eo,i) = 
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Function F{X) in formula (1) (scaling function) is normal- 
ized LDF with respect to variable X = r/R^,s,. It does 
not depend practically on primary particle type, shower en- 
ergy and age. Besides it is not sensitive to variations of basic 
parameters of hadronic interaction model. Finally we ob- 
tain the following formula for electron density (A. A. Lagu- 
tin, R.l. Raikin, 2001): 
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Though the structure of function (3) is similar to modified 
Linsley function, which is traditionally used for fitting ex- 
perimental LDF of charged particles measured by scintilla- 
tion counters, the important difference is that function (3) is 
one-parametric and scale-invariant. 

Our latest results obtained for the mean square radius of 
electrons can be approximated as follows (A. A. Lagutin et 
al., 2001; R. 1. Raikin et al., this proceedings): 
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Here p{t) is air density at depth t, po = 1.225 g/cm^, 
^max - depth of maximum of average cascade curve. Ap- 
proximation (4) gives one- valued relation between f^ax and 
normalized average lateral distribution function of electrons. 
It is important that relation (4) is stable to variation of 
hadron-air inelastic cross section and inclusive spectra of 
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Fig. 1. Average lateral distribution of electrons for 10 eV proton- 
induced extensive air showers. Function (3) is shown in comparison 
with CORSIKA/QGSJET data (M. Nagano et al., 2000) 



secondaries remaining valid within the limits defined by 
widespread hadronic interaction models (A. A. Lagutin et al., 
2001; R. I. Raikin et al., this proceedings). 

On Fig. 1 we compare our scaling LDF of electrons with 
CORSIKA/QGSJET predictions (M. Nagano et al., 2000) for 
10^^ eV proton-initiated extensive air showers. Our results 
are presented by function (3) with = 4.37 • 10* and 
i?m.s. = 118m(Uax = 722 g/cm2) as QGSJET predicts for 
proton-induced showers of considered energy. It is seen, that 
function (3) is steeper than CORSIKA data. The difference 
becomes essential at r > 1 km and increases considerably 
when irregularities appear in CORSIKA data due to the lack 
of distant particles caused by thin samphng. 



3 Comparison of electron LDF witli AGASA data 

Adequate comparisons of theoretical results with experimen- 
tal lateral distributions of charged particles require to take 
into account the contribution of low-energy muons and also 
detector response effect. 

According to CORSIKA simulations (M.Nagano et al., 
2000) lateral distribution of muons changes its form sharply, 
when muon threshold energy decrease from 1 GeV to 
0.25 GeV. About 10-15% additional increasing in muon den- 
sities without disturbing the shape of LDF was pointed out 
for threshold energy 10 MeV. As a result coefficient fc^(r) = 
Pij.i'f": > 10 M;eV)/p^(r; > 1 GeV) in vertical showers in- 
creases from ~ 1 . 8 at 600 m to ~ 2 . 7 at 2000 m from the core 
making muon lateral distribution much flatter in comparison 
with LDF measured experimentally by muon detectors with 
high threshold energy. 

In several recent papers (T. Kutter, 1998; A. A. Lagutin 
et al., 1999; A. 1. Goncharov et al., 2000; M.Nagano et al.. 
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2000) theoretical study of the influence of scintillation detec- 
tor response on the shape of charged particle LDF has been 
made. The correction factor ksc{r) = S{r)/pch{r), where 
S{r) is scintillation yield in units of minimum ionizing par- 
ticle {scintillator density) and Pch{r) is charged particle den- 
sity itself, was investigated. While absolute values of kgc 
calculated using different methods and simulation codes are 
different from each other, the main tendency in radial de- 
pendence is same. Function ksc{r) increases distinctly from 
several tens to at least several hundreds meters from the core 
position. It means that utiUzing a factor 1.1 of scintillator 
density to spark chamber density, which has been determined 
experimentally by AGAS A group within 100 m from the core 
(M. Teshima et al., 1986), is not quite accurate for large dis- 
tances. For our comparisons we implemented data T. Kutter 
(1998), calculated by CORSIKA for 5 cm plastic scintillator, 
according to which kgc ~ 1.1 at r = 50 m and ksc ~ 1-4 
at r = 600 — 1000 m. For very large distances (r > 1 km) 
no changes in the shape of LDF due to scintillation detector 
response were assumed. 

Formally, the relation between two experimentally observ- 
able characteristics S{r) and /3^(r; > 1 GeV) can be pre- 
sented as follows: 

^ [Pe(r) + p^(r;>lGeV)fc^(r)]fc,e(r) _ 

-^array 

Here -f^array IS Coefficient, which depends on single particle 
definition used in concrete experiment [in case of AGASA 
K^„a.y = 1-1 (M.Nagano et al., 2000)]. 

Since k^{r) and ksc{r) are determined basically by low 
energetic part of a cascade, it is reasonable to expect that ra- 
dial behaviour of both correction factors will not depend crit- 
ically on features of hadronic interactions at very high ener- 
gies, primary energy (at least in a limited energy range) and 



composition. This conclusion has been partially confirmed 
by theoretical studies T. Kutter (1998); A. I. Goncharov et al. 
(2000); M. Nagano et al. (2000). Therefore we can use rela- 
tion (5) to estimate electron LDF from the experimental data. 
On Fig. 2 we show the result of such estimation of pdr) 
for vertical 10^* eV showers obtained from AGASA data 
in comparison with our scahng distribution for proton pri- 
maries (same as on Fig. 1, except of normalizing constant). 
AGASA empirical functions for S{r) and /9p(r; > 1 GeV) 
together with estimated p^{r; > 10 MeV) are also shown. 
We utilized both logarithmic (a) and hnear (b) radial distance 
scales to emphasize that our result is in very good agreement 
with the one estimated from AGASA experimental data in 
whole radial distance range. 

4 Cosmic ray primary composition from the shape of 
LDF 

The chemical composition of primary cosmic rays with Eq > 
10^^ eV is not measured well for today. The basic method 
used to deduce primary composition from experimental data 
of large ground-based air shower arrays is muon compo- 
nent analysis. Though muon local density at fixed distance 
from the core position is rather sensitive to the type of pri- 
mary particle (including 7-rays), the observable muons with 
Eth = (0.5 — 1) GeV cos are mostly derived from high 
energy interactions. Thus the interpretation of experimen- 
tal data on the basis of comparisons with theoretical results 
depends on hadronic interaction model used in calculations. 
This problem is still unsolved. 

A key result, which allows us to consider the shape of 
lateral distribution function measured by scintillation coun- 
ters as a source of information about primary composition 
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Fig. 3. The energy dependence of the mean square radius of EAS 
electrons for proton- and iron- initiated showers at AGASA obser- 
vation depth as predicted by approximation (4) with CORSIKA and 
SIBYLL results for tmax. See text for details 

of ultrahigh energy cosmic rays, is rather strong energy (and 
consequently composition) dependence of the shape of elec- 
tron LDF. It is important that, as scaling property (1) shows, 
this dependence can be completely described by the variation 
of single parameter - mean square radius of electrons. Ac- 
cording to our results (A. A. Lagutin et al., 2001; R.I. Rai- 
kin et al., this proceedings) the narrowing of electron LDF 
with energy is likely to be conserved at least up to 10^° eV, 
while the dependence on basic parameters of hadronic inter- 
action model is relatively weak. On Fig. 3 we show mean 
square radius of electrons vs. primary energy for proton- and 
iron-initiated showers at AGASA observation level (tobs = 
920 g/cm^). The shaded areas represent uncertainties of the 
result among QGSJET (top) and SIBYLL (bottom) interac- 
tion models [here we used approximation (4) with these mod- 
els predictions for t,nax as reported in A. M. Hillas (1997)]. 

Unfortunately, the mean square radius itself can not be 
estimated directly from experimental data so far as it de- 
mands high precision measurements of lateral distribution 
of electron component in wide radial distance range (from 
several meters to several thousands meters from the core). 
Nevertheless, the change of slope of charged particle LDF 
around 600 m from the shower core reflects well the changes 
in electron lateral distribution, as experimental uncertainties 
in this region are relatively small and the contribution of 
muons with Eq^ — 1 GeV in vertical showers is about 10% 
at 5 ■ 10^^ eV, decreasing with energy (N.Hayashida et al., 
1995). 

We should point out that theoretical studies of both fc^(r) 
and ksc{r) are far from being complete up to now and also 
it is still difficult to finally exclude some exotic processes or 
new particles in EAS of ultrahigh energies. Furthermore, it 
is impossible to confirm or deny the presence of small frac- 
tion of 7-induced showers from such kind of analysis. But 
the main disadvantage for making some conclusions about 
composition is controversial character of existing experimen- 



tal data, obtained by AGASA and Yakutsk array. While en- 
ergy independence of the shape of S{r) from ^ 10^^ ''' eV up 
to highest observed energies reported by AGASA group re- 
peatedly [see, for example, N. Hayashida et al. (1999)] is an 
evidence for changing primary composition from relatively 
tight to relatively heavy in this energy region, Yakutsk data 
shows extremely sharp steepening of LDF around 5 • 10^* eV 
(A. V. Glushkov et al., 1997, 1999), which (refusing extraor- 
dinary explanations) leads to converse conclusion. Neverthe- 
less the weak sensitivity of the result to hadronic interaction 
model makes such kind of analysis interesting as another one 
source of information about primary composition, hopefully 
in nearest future with further accumulation of experimental 
data. 
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